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We investigate the origin of order in the low-lying spectra 
of many-body systems with random two-body interactions. 
Contrary to the common belief our study based both on an- 
alytical as well as on numerical arguments shows that these 
are the higher J-sectors whose ground states are more orderly 
than the ones in the J = 0-sector. A predominance of J = 
ground states turns out to be the result of putting on together 
states with different characteristic energy scales from different 
J-sectors. 

PACS numbers: 05.45. -l-b, 21.60.Cz, 24.60.Lz 

The nature of mechanism generating order out of ran- 
domness constitutes one of the most fundamental issues 
of the contemporary physics. Theories based on vari- 
ous versions of ensembles of the random matrices pro- 
vide one possible theoretical frame for studying such ef- 
fects. In this context the recently identified prepon- 
derance of the J = ground states in strongly interacting 
Fermi systems, such as atomic nuclei, arising from ran- 
dom two-body interactions came as a surprise since there 
is no obvious pairing character in the assumed random 
force. Various possible explanations of this effect have 
been tested ^ with no success, however. One purpose 
of the present note to provide a consistent picture of the 
mechanism generating this effect. 

Schematically, indicating nevertheless all the relevant 
ingredients, the interaction matrix elements ^, of good 
total angular momentum J in the shell-model basis \a) 
can be expressed as follows B: 
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The summation runs over all combinations of the two- 
particle states \i) coupled to the angular momentum J' 
and connected by the two-body interaction g. g^^, denote 
the radial parts of the corresponding two-body matrix 
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globally represent elements of the 



elements while 

angular momentum recoupling geometry. 

In statistical ensembles of matrices the crucial factor 
determining the structure of eigenspectrum is the prob- 
ability distribution Pv(v) of matrix elements |^]. Es- 
pecially relevant are the tails of such distributions since 
they prescribe the probability of appearance of the large 
matrix elements. From the point of view of the mecha- 
nism producing the energy gaps they are most effective 
in generating a local reduction of dimensionality respon- 
sible for such effects. In principle, the probability distri- 
bution of the shell model matrix elements is prescribed 
by their general structure expressed by the eq. (0), pro- 



are known. In general terms this structure can be con- 
sidered to have the following form: 



(2) 



and each Vi to be a product of another two variables 
denoted as Ci and Gi. By making use of the convolution 
theorem Q the probability distribution Pv{v) that V 
assumes a value equal to v can be expressed as: 



Pv{v) = F-^[F{PvM)) ■ FiPvM) 



FiPvJvN))], 
(3) 



where F denotes a Fourier transform, F^^ its inverse and 
Pvii^i) the probability distributions of individual terms 
in eq. (S). Taking in addition into account the fact that 



PvM)= I dg,PG^{g,)PcA-)^, 
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one can explicitely derive the form of Py(w) in sev- 
eral cases. Assuming for instance that all the above 
constituents are identically Gaussian distributed (then, 
according to eq. (P, P v-(fi) ~ KQ{\vi\)/'K and thus 
F{Pvi{vi)) = l/vT+tJ^ ) one arrives at 



Pv{v) 
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where K stands for the modified Bessel function. Asymp- 
totically, for large v, this leads to 



Pviv) ^ exp{-\v\) \vf/^-\ 
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vided the probability distributions of both g'/^, and c 
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For such a global estimate the identical Gaussian dis- 
tribution of g^^, is consistent both with the Two-Body 
Random Ensemble (TBRE) Q and with the Random 
Quasiparticle Ensemble (RQE) The only anticipated 
difference originates from the fact that in the second case 
the variance of the distribution drops down with J' like 
the inverse of 2 J' -I- 1 which is expected to result in a 
smaller effective N as compared to TBRE. By contrast, in 
both versions of the above random ensembles the geom- 
etry expressed by Cj"" enters explicitely. However, the 
complicated quasi-random coupling of individual spins is 
believed [Q to result in the so-called geometric chaotic- 
ity ^ . If applicable indeed then this fact implies validity 
of the above estimate for Pv{v). Since this is an im- 
portant element for the present considerations below we 
explicitely verify its range of applicability. 

The model to be quantitatively explored here consists, 
similarly as in ref. Q], of 6 identical particles (all single- 
particle energies are set to zero) operating in the sd shell. 
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From the nuclear spectroscopy point of view this can be 
identified as ^^O. Fig. 1 shows distributions of the cor- 
responding geometrical factors Cj"" for a ^ a' and for 
the relevant values of J and J'. 
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Fig. 1 The normalised distribution of geometrical factors 
Cj"" entering the off-diagonal matrix elements (eq. (|^)) 
for the model of 6 particles in the sd-shell. 



As one can see, the Gaussian may be considered quite a 
reasonable representation of the distribution of such fac- 
tors for all combinations of J and J' with one exception 
for those which involve J = 0. In this later case the dis- 
tribution of Cj"" resembles more a uniform distribution 
over a finite interval located symmetrically with respect 
to zero. These empirical facts justify well the estimates 
of Pv{v) based on eq. (j^) for J ^ Q and not so well for 
J = 0. More appropriate in this particular case is to 
assume a uniform distribution of Cj?j" over an interval 
confined by say — co and co, i.e., Pci(ci) = l/2co, retain- 
ing PQ.[gi) in its original Gaussian form of course. By 
making use of eqs. (0) and (0) one then obtains 
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which for large v behaves like 
Pv{v) ~ exp(- 



(8) 



An explicit calculation of the distribution of the shell 
model off-diagonal matrix elements for the various J- 
values based on the present model with two-body matrix 
elements drawn from RQE (TBRE results in similar rela- 
tions among different J-sectors though the distributions 
are somewhat broader as compared to RQE) confirms the 
above analytical estimates as is illustrated in Fig. 2. 
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Fig. 2 The probability distribution of nonzero many- 
body off-diagonal matrix elements in different J-sectors 
drawn from one thousand of RQE samples of two-body 
matrix elements. The energy scale is set by v where 
wj' = /{2J' + 1) and wji determines the RQE mean 
square variance. These distributions are fitted (continu- 
ous solid lines) in terms of eq. with N treated as a 
fitting parameter. The corresponding best iV's (Neff) for 
each J are listed. By increasing N the distribution pre- 
scribed by eq. (H) quickly approaches (as a consequence 
of the central limit theorem) the Gaussian distribution. 
In this way the J = distribution is demonstrated to 
be much closer to the Gaussian than the remaining ones 
whose asymptotic behavior is consistent with a slower, 
exponential fall-off. 
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Indeed, such a distribution in the J = sector resem- 
bles more a Gaussian and the tails of this distribution 
drop down faster as compared to the J ^ 0-sectors where 
the large v tails drop down slower, as consistent with an 
exponential asymptotics of eq. (^). At the same time 
the J 7^ sectors are dominated by very small matrix 
elements to a larger degree than J = 0. The probabil- 
ity of appearance of a large off-diagonal matrix element 
which in magnitude overwhelms the remaining ones is 
thus greater for J 7^ than for J = 0. Such an effective 
reduction of the rank in the former case is expected to 
result in a stronger tendency to localization as compared 
to GOE The corresponding characteristics can be 

quantified in terms of the information entropy 

a=l 

of an eigenstate labelled by I from the J-sector. The 
coefficients a{ ^ denote the eigenvector components in the 
basis |a|. Since the definition of Kf involves the total 
number of states Mj which differ for different J's and 
since the appropriate reference for our considerations is 
GOE, we normalise Kf to the GOE limit of this quantity 
which reads pX| ]: 

KhoE = ^{M.j/2 + 1) - VX3/2), (10) 

where ijj is the digamma function. Within our model 
the so-calculated and RQE ensemble averaged quantity 
for all the states versus their corresponding energies Ei 
is illustrated in Fig. 3. As anticipated, it is not J = 
whose lowest eigenstate comes out most localised, i.e., 
most regular. The lowest states for several higher J val- 
ues deviate much more from GOE. This in particular 
applies to J = 2 and, especially, to J = 4. This thus 
indicates more favorable conditions for the emergence of 
energy gaps for larger J than for J = 0, contrary to the 
numerical outcome of ref. ||l[ . That our numerical proce- 
dure is equivalent to that of ref. is confirmed by the 
fact that we can reproduce the results when calculating 
the same quantities. The question thus arises as how to 
reconcile such conflicting conclusions. 

In fact Fig. 3 provides one hint of what might be the 
reason. The J = states are spread over the broadest 
energy interval even though the number of states (Mq = 
14) is here significantly smaller than for several larger 
J values (Ml = 19, M2 = 33, M3 = 29, M4 = 26). As 
a result, the average separation between the states is a 
factor of few larger for J — Q than for the remaining ones. 
Putting on the states from different J-sectors together, 
as is done in ref. Q], is thus likely to hide the genuine 
character of the effect under consideration. 

In order to illustrate the consequence of such a proce- 
dure in Fig. 4 (dashed line) we show distributions of the 
ground state [E() gaps 

= {Ei-E()/D-\ (11) 
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Fig. 3 The information entropy normalised to its GOE 
limit {Kj' / Kqq^) for all the states I from various J- 
sectors (all positive parity) versus energies (-E/) of those 
states. All the quantities are ensemble averaged. The 
energy units are the same as in Fig. 2. 

where similarly as in ref. [Q| 

D' =<Ei-Ei>, (12) 

for each J individually. 

As it is clearly seen the J — 0-sector no longer distin- 
guishes significantly from the remaining ones. From the 
point of view of our considerations presented above one 
would however expect an even reduced probability for 
occurance of the large ground state energy gaps in this 
particular sector as compared to the J > ones. As the 
solid lines in Fig. 4 indicate such an effect does indeed 
take place when in eq. (Ill]) is replaced by 

D' =< Ei,^ - Ei > /{Mj - 2). (13) 

Actually the definition of D'^ as specified by eq. (|l2|) looks 
somewhat arbitrary. It seems more appropriate and more 
consistent with the above global considerations to relate 
the ground state energy gap just to the average global 
level spacing among the remaining states, characteristic 
for a given J. 

Finally one may ask a question why this tendency does 
not extend to even higher J-values. In this connection 
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one has to remember that the ofF-diagonal matrix ele- 
ments is not the only relevant element. These are the 
diagonal matrix elements which constitute the driving 
term. Irrespective of the value of J their distribution is 
always Gaussian-like. This can be observed numerically 
and is consistent with arguments formulated in terms of 
eqs. (0 - since the geometrical factors Cj"" entering 
the diagonal matrix elements are always nonnegative. In- 
creasing however J beyond 4 results in a significant re- 
duction of the variance of Pv{v) for the diagonal matrix 
elements and consequently a larger fraction of the off- 
diagonal matrix elements becomes effective in mixing the 
basis states. Thus the effect under investigation results 
from an interplay between distributions of the diagonal 
and off-diagonal matrix elements and their relative mag- 
nitudes. 




Fig. 4 Distributions of ground state energy gaps as 
defined by the eq. (|ll|) for successive J's. The dashed 
line uses defined by eq. while the solid line the 
one defined by eq. (p^). 

In conclusion, the present investigation based both on 
theoretical as well as on numerical arguments clearly 



shows that the many-body problems described in terms 
of various variants of the two-body random ensembles 
(like RQE or TBRE) develop quantitatively well identi- 
fied deviations from the GOE. These deviations quan- 
tified in terms of localisation or of energy gaps point 
to the angular momenta between 2 and 4 as those J- 
sectors whose ground states are ordered most. From 
this perspective the current interpretation of the numer- 
ical results of ref. |^] is thus an artifact of mixing states 
with different characteristic energy scales from different 
J-sectors. As a side remark it seems also appropriate to 
notice at this point that the arguments going in parallel 
to eqs. provide a more adequate approach towards 

understanding the distribution of matrix elements in real- 
istic nuclear shell-model calculations than the ones based 
on multipole expansion 
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